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Strongly correlated analogues of topological insulators have been explored in systems with purely on-site 
symmetries, such as time-reversal or charge conservation. Here, we use recently developed tensor network tools 
to study a quantum state of interacting bosons which is featureless in the bulk, but distinguished from an atomic 
insulator in that it exhibits entanglement which is protected by its spatial symmetries. These properties are 
encoded in a model many-body wavefunction that describes a fully symmetric insulator of bosons on the honey¬ 
comb lattice at half filling per site. While the resulting integer unit cell filling allows the state to bypass ‘no-go’ 
theorems that trigger fractionalization at fractional filling, it nevertheless has nontrivial entanglement, protected 
by symmetry. We demonstrate this by computing the boundary entanglement spectra, finding a gapless entan¬ 
glement edge described by a conformal field theory as well as degeneracies protected by the non-trivial action 
of combined charge-conservation and spatial symmetries on the edge. Here, the tight-binding representation of 
the space group symmetries plays a particular role in allowing certain entanglement cuts that are not allowed on 
other lattices of the same symmetry, suggesting that the lattice representation can serve as an additional symme¬ 
try ingredient in protecting an interacting topological phase. Our results extend to a related insulating state of 
electrons, with short-ranged entanglement and no band insulator analogue. 


I. INTRODUCTION 

The experimental discovery [1] of topological band insu¬ 
lators that cannot be adiabatically continued to the atomic 
limit as long as time-reversal invariance is preserved [2, 3] 
has spurred the exploration of a broad array of phases where 
symmetries protect subtle, non-local features that distinguish 
them from trivial, unentangled insulators. These phases, 
collectively known as symmetry-protected topological (SPT) 
phases [4], have by now been observed in several experimen¬ 
tal realizations in one, two and three dimensions and an ex¬ 
tensive mathematical framework has been developed for their 
characterization and classification [5-11]. 

Although a system that is in such a phase is often only 
subtly distinguished from a trivial insulator in its bulk, the 
topological properties can be observed by focusing on the 
boundary of a finite system. In many cases, the boundary 
exhibits gapless features such as the localized Majorana zero 
modes at the ends of one-dimensional topological supercon¬ 
ductors [12], the helical edge states of the quantum spin Hall 
effect [13-15] and the protected Dirac cones on the surface of 
three-dimensional Z 2 topological insulators [16]. It has been 
shown that these features are also manifest in the entangle¬ 
ment spectrum [17], where they generally take the form of 
either protected degeneracies or gapless spectra that mirror 
the gapless modes on a boundary. These signify entanglement 
that cannot be removed while preserving the symmetry, and 
can thus be used to establish that SPT phases are fundamen¬ 
tally distinct from trivial, non-entangled insulators. 

Considering not only local symmetries but also symmetries 
that relate the physical locations of degrees of freedom, such 
as the spatial symmetries of rotation or reflection, can lead to 
an even richer panoply of phases [18-23] that require new 
probes as well as a more powerful classification. The role 


of the boundary is also modified when spatial symmetries are 
considered; while in some cases, any possible physical bound¬ 
ary of the system may break the relevant symmetry, the non¬ 
trivial features can still be extracted from the entanglement 
spectrum. This allows the (numerical) identification of such 
SPT phases even when a necessary ingredient for the topo¬ 
logical protection is a non-local symmetry operation. How¬ 
ever, while non-interacting topological crystalline insulators, 
i.e. band insulators whose non-trivial structure is protected by 
crystal symmetries, have been well explored [24, 25], far less 
is known about their interacting counterparts in two or three 
dimensions. 

In this work, we compute the entanglement properties of an 
insulating state of interacting bosons on the honeycomb lat¬ 
tice, and show that it constitutes a topological phase protected 
by lattice symmetries. In particular, we show that the non¬ 
trivial entanglement is not related just to the group formed by 
the lattice symmetries, but becomes tied to the specific real¬ 
ization as a honeycomb lattice. Combined with the symmetry, 
this enforces a non-trivial short-range entanglement structure. 

The wavefunction we consider is an insulator of interact¬ 
ing bosons on the honeycomb lattice at a Ailing of one boson 
per unit cell, or half Ailing per site. It forms one example of 
a class of insulators that require a non-Bravais lattice, i.e. a 
lattice with multiple sites per unit cell. The necessity for such 
a non-trivial unit cell arises due to higher-dimensional gen¬ 
eralizations [26-29] of the Lieb-Schultz-Matthis (LSM) the¬ 
orem [30],which forbids the existence of a featureless state 
— a state that neither spontaneously breaks a symmetry, nor 
displays intrinsic topological order, nor has power-law corre¬ 
lations and is thus “gapless” — in systems with a fractional 
Ailing per unit cell. While such a featureless state at half- 
fllling per site is allowed on the honeycomb lattice, the explicit 
construction is challenging. Symmetry guarantees that a free- 
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fermion spectrum is gapless at certain high-symmetry points, 
and there is thus no basis of localized, symmetric and orthogo¬ 
nal Wannier states. This implies that a featureless state on the 
honeycomb lattice cannot be constructed by filling a perma¬ 
nent of localized Wannier orbitals [31], and any construction 
of a quantum state thus must involve interactions. Ref. 32 pur¬ 
sued an approach of constructing a permanent wavefunction 
by filling local and symmetric orbitals that are not orthogonal 
and it was argued that that the resulting wavefunction is in¬ 
deed featureless. In particular, using numerical simulations it 
was found that the state exhibits isotropic and exponentially 
decaying correlations, and arguments were presented that it is 
not topologically ordered. 

Here while we confirm the featureless bulk of the state, in¬ 
cluding the absence of intrinsic topological order, we show 
that nevertheless the entanglement of the state cannot be re¬ 
moved while preserving certain symmetries — it is symmetry- 
protected. The relevant symmetry is a combination of charge 
conservation and lattice symmetries, which together protect 
universal features in the entanglement spectrum. In particu¬ 
lar, we show that the low-lying entanglement spectrum is to 
great accuracy described by that of a c = 1 conformal field 
theory, and that there is an exact double degeneracy through¬ 
out the entanglement spectrum for certain geometries, which 
is protected by the symmetries of the state and thus serves as 
a topological invariant identifying the SPT order. Since lat¬ 
tice symmetry is involved crucially, this provides one of the 
first examples for an SPT of interacting bosons protected by 
lattice rather than on-site symmetries. To further substantiate 
the robustness of these features, we obtain parent Hamilto¬ 
nians for the phase in certain quasi-one-dimensional geome¬ 
tries and study the effect of weak symmetry-preserving and 
-breaking perturbations on the ground states of these Hamil¬ 
tonians. 

All of these features become accessible through a de¬ 
scription of the state as a projected entangled-pair state 
(PEPS) [33-36]. These states form a specific class of ten¬ 
sor network states that corresponds to a generalization of the 
well-known matrix-product state (MPS) [37^0] framework 
to higher dimensions. PEPS have been shown to be a power¬ 
ful description of many classes of gapped systems, including 
topologically ordered and SPT phases. Here, we have an exact 
description of the state as a PEPS, allowing us to extract prop¬ 
erties such as the entanglement spectrum and the topological 
invariants exactly on certain geometries; we emphasize that 
these properties of the state are not accessible to other numer¬ 
ical methods. 

The topological invariants extracted here form examples 
of a broad class of invariants that provably must be constant 
throughout the phase. These differ from the order parameters 
that measure local symmetry breaking in that they are not re¬ 
lated to the expectation values of local operators. Early exam¬ 
ples of topological invariants for SPT phases are the string or¬ 
der parameter for the one-dimensional AKLT phase [22, 41- 
43], and the spin Chern number for the quantum spin Hall 
effect [44]. The invariants we consider here measure how the 
action of the symmetry is implemented on the physical edge 
states of open systems or on the Schmidt states of an entan¬ 


glement decomposition [45]. These invariants feature heavily 
in the classification of SPT phases with on-site symmetry, and 
similar invariants that apply to free-fermion states have been 
used for topological crystalline insulators [24, 46]. In con¬ 
trast, topological invariants for interacting states protected by 
lattice symmetries in more than one dimension are poorly un¬ 
derstood. We will discuss the action of the symmetry on the 
edge of the state and progress towards the goal of finding a 
topological invariant to identify the corresponding phase. 

The rest of this paper is structured as follows: in Sec¬ 
tion II, we review the honeycomb featureless boson insula¬ 
tor (HEBI) state, and introduce its PEPS representation. In 
Section III, we discuss results for the correlation functions of 
this state. In Section IV, we discuss the entanglement spec¬ 
tra that we obtain numerically from the PEPS representation 
and discuss in detail their connection to the spectrum of a 
free boson conformal field theory. In Section V, we describe 
the symmetry-protected topological invariants that allow us 
to identify the symmetries that protect certain entanglement 
properties of the state in quasi-one-dimensional geometries. 
Section VI discusses the effect of weak perturbations to a par¬ 
ent Hamiltonian in the same quasi-one-dimensional setup, and 
Section VII introduces a different perspective on the phase 
from the point of view of a boson-vortex duality. 

II. CONSTRUCTION OF THE FEATURELESS BOSON 
INSULATOR 

In the honeycomb lattice, each unit cell is associated with 
exactly one hexagon plaquette, which respects the lattice point 
group symmetries. As shown in Ref. 32, this provides an ex¬ 
plicit construction of a bosonic insulator on the honeycomb 
lattice that is completely featureless in the bulk, henceforth 
referred to as honeycomb featureless boson insulator (HEBI). 
The state is succinctly described by the following expression: 

o ViGO / 

Here, O denotes the elementary hexagons of the honeycomb 
lattice. Despite the deceivingly compact expression, this 
many-body bosonic state is strongly correlated and requires 
concrete computation for its properties to be unveiled. 

We focus on two closely related variants of this state: a ver¬ 
sion of soft-core bosons where b\ creates a boson on site i and 
obeys the usual bosonic commutation relations, and a hard¬ 
core version of the same state where bj also creates a boson 
but {bl)'^ = 0. In either case, the operator creates 

exactly one boson per hexagon; as there is one hexagon per 
unit cell of two sites of the lattice, the state has one boson 
per unit cell, or half a boson per site, thus allowing the exis¬ 
tence of a featureless state. In the case of soft-core bosons, the 
maximum number of bosons per site is three. 

Ref. 32 examined properties of both the soft-core and hard¬ 
core variants of this state. In the soft-core case, ground-state 
correlations were mapped to those of a classical loop model on 
the triangular lattice. A Monte Carlo analysis thereof revealed 
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Figure 1. Floneycomb lattice PEPS and zig-zag entanglement cut. 
In this PEPS of rank-4 tensors, the top and bottom edges are iden¬ 
tified, forming a cylinder with circumference FK = 3 unit cells. A 
one-dimensional MPS representation is constructed by contracting 
the tensors in each cylinder slice (region marked by dotted lines). 
The entanglement cut used (either one of the dotted lines) passes 
through the hexagon mid-points, preventing the tight-binding lattice 
from gaining additional sites as long as crystalline symmetries are 
preserved. 

that the boson Green’s function {b\b^) decays exponentially 
- thereby ruling out the possibility that the many-body wave- 
function describes a superfluid - and further, that the corre¬ 
lation functions of a variety of neutral operators (e.g., those 
describing charge- or bond-density order) remained short- 
ranged. This loop model mapping also included a variational 
parameter, m, that tunes the the soft-core boson wavefunction 
on the honeycomb lattice into that of a trivial Mott insulator 
on a triangular lattice of fictitious sites placed at the center of 
each hexagon. The absence of a transition under this perturba¬ 
tion was taken as evidence that the ground state would remain 
unique on manifolds of nontrivial topology, thereby ruling out 
the possibility that the wavefunction describes a topologically 
ordered phase. For the case of hard-core bosons, a different 
quantum-classical mapping enabled the efficient calculation 
of boson number correlations. This directly revealed that the 
hard-core projection did not induce any long-range correla¬ 
tions in the neutral sector. Although working in the num¬ 
ber basis precludes direct access to ‘charged’ correlators such 
as the boson Green’s function, on general grounds, the alge¬ 
braic decay characteristic of classical 2D superfluids [47] is 
expected to also infect density-density correlations, and thus 
their exponential decay provides indirect evidence that the 
hard-core boson wavefunction also lacks superfluid order. We 
note that none of these quantum-classical mappings can read¬ 
ily provide insights into the entanglement properties of the 
wavefunction. 


A. PEPS representation 

In order to go beyond the properties accessible via these 
quantum-classical mappings of the HFBI, and in particular in 
order to be able to study its edge properties, we now derive a 
representation as a projected entangled pair state (PEPS). Im¬ 
portantly, this PEPS description will respect all of the relevant 
symmetries of Eq. (1). 


Figure 2. Intermediate tensor network for HFBI state. Here, the 
tensors labeled D are located on the sites of the honeycomb lat¬ 
tice, while the tensors labeled W are located on the centers of each 
hexagon. Dotted lines thus represent the physical lattice, while the 
solid lines indicate auxiliary bonds over which the tensor network is 
contracted. 

To obtain a PEPS construction, we first choose a local basis 
|n) of boson occupation numbers, i.e. b%\n) = n\n). The 
PEPS will thus describe the coefficients of |'0) in this basis, 
(ni... nil'!/)). The PEPS representation is most easily ob¬ 
tained in a two-step construction, where we first construct the 
state shown in Pig. 2. Here, the tensor labeled W = 
which is placed in the center of each hexagon, is a rank-6 ten¬ 
sor given by 

f 1 : = 1 

^ (2) 

0 : else 

where each Ui G {0,1}. 

This tensor describes the coefficients of a so-called W- 
state in the occupation number basis, i.e. _ 

(ni .. .ne \ 10)- We note that this tensor is symmet¬ 

ric under permutations of its indices. 

On the sites of the physical lattice, we have placed a rank-4 
tensor denoted as D, shown in panel (a) of Pig. 3(a), which 
connects the W tensors from three adjacent hexagons, and as 
fourth index has a physical index p. Por a state of soft-core 
bosons, where p = 0,1, 2, 3, this tensor is given by 

^sc ^ f VW- ■ P = io + ii + i2 .os 

I 0 : else ’ ^ ’ 

We can also encode a state of hard-core bosons by replacing 
Dhy 

T^hc ^ / 1 : P = *0 + H + *2 < 1 ..s 

I 0 : else ' ^ ’ 

Other values for the D and W tensors that respect the charge 
and lattice symmetries can also give rise to featureless insula¬ 
tors. Some of these variants are described in Appendix C. 

This tensor network wavefunction manifestly respects all 
the translational and point group symmetries of the honey¬ 
comb lattice, since the tensors W and D are invariant under 
rotations of their virtual indices in the plane. One can also 
check that the wavefunction is (7(1) invariant with charge 1 
per plaquette. 
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(a) D tensor 


= 

(b) W-tensor and factored form 



lattice, it is an exact representation of the FBI state and does 
not break any symmetries after contracting the indices. 

This decomposition respects the physical 17(1) charge con¬ 
servation symmetry in that all tensors are separately C/(l)- 
invariant [48]. To make this manifest, we have indicated in 
Fig. 3 arrows on each bond that show the flow of charge. 




(c) A tensor (d) B tensor 



(e) PEPS tensor network for F.B.I. state 

Figure 3. Construction of PEPS for HFBI state. The site tensors 
(shown in panels (c), (d)) are constructed using the factors of the pla- 
quette tensor W (panel (b)) combined with the original vertex tensor 
D (panel (a)). The red line in panel (e) shows where the entangle¬ 
ment cut considered in this paper crosses the network. 


In order to convert the tensor network of Fig. 2 into a PEPS 
representation, we first factor the FF-tensor into a matrix- 
product state of six tensors as shown in Fig. 3(b). We choose a 
form of the MPS that breaks the rotational symmetry of the W- 
state (which appears as translational symmetry of the MPS). 
This allows us to obtain an MPS description with a small bond 
dimension of M = 2; a fully symmetric choice would require 
bond dimension 6. Since these states are physically equiva¬ 
lent, all physical quantities are unaffected by this choice. One 
possible decomposition is given by 


W 




= E 


T/*1 FF*2 FF*5 X*® 

^ CL\' ^ a\Ci'2 04^5 Q5 


( 5 ) 


where and 


, :*o+i = *i 

1 0 : else 


where each index takes values in {0,1}. Applying this to each 
FF-tensor yields the state as shown in Fig. 3(e). By contract¬ 
ing the four tensors in each shaded region together, we obtain 
a PEPS in the regular form as shown in Pig. 1. The result¬ 
ing PEPS has a bond dimension of M = 2 on the horizontal 
bonds, and a bond dimension of M = 4 on all other bonds. 
While it superficially breaks the rotational symmetry of the 


B. Representation on infinite cylinders 

Por the calculations presented in this manuscript, we con¬ 
sider the state |^) on a cylinder of infinite length, but finite 
circumference of FF unit cells. In Fig. 1, we have indicated 
the choice of boundary conditions for the cylinder used in this 
paper. For many practical purposes, the PEPS on an infinite 
cylinder can be represented as an infinite, translationally in¬ 
variant matrix-product state of bond dimension x = 2'^ and 
physical dimension p = 4^^ (p — 2^^) for the soft-core 
(hard-core) case. The MPS is created by blocking all tensors 
in each slice of the cylinder, as shown in Figure 1. 

With each cylinder slice blocked together and considered 
as an MPS, the procedures we use for computing both corre¬ 
lation functions and entanglement properties are in principle 
identical to those used previously in MPS [45, 49]. Due to the 
exponential increase in the MPS bond dimension, this numeri¬ 
cally exact approach scales exponentially in the circumference 
of the cylinder. It is however computationally advantageous 
to exploit the additional structure present in the PEPS transfer 
operator; by doing so, we can compute correlation functions 
and the entanglement spectrum for the cut shown in Figure 
3(e) for the HFBI state on cylinders of circumference up to 
FF = 10. These computations are presented in the following 
Sections III and IV. 


III. CORRELATION FUNCTIONS 

In Ref. 32, certain real-space correlation functions of the 
featureless boson insulator state were studied using a mapping 
to a particular classical statistical mechanical system which 
was sampled using Monte Carlo techniques. Here, we go be¬ 
yond this by employing PEPS calculations on infinite cylin¬ 
ders that allow us to measure a broader class of correlation 
functions and, in particular, allow us establish a strict upper 
bound on the exponential decay of al/ two-point correlation 
functions for an infinite cylinder of given width. 

In Fig. 4, we show both density-density and off-diagonal 
short-range correlation functions for a cylinder of circumfer¬ 
ence FF = 8. Comparing these to the Monte Carlo results 
of Ref. 32, which have been computed for a different geome¬ 
try, we And good qualitative agreement. Crucially, while the 
boundary conditions we choose break the rotational symmetry 
by making the system periodic in one direction and infinite in 
the other, the short-range correlations for distances up to half 
of the cylinder circumference appear unaffected by this. 

It is a well-known result that PEPS can, in the thermo¬ 
dynamic limit, exhibit power-law correlation functions [50], 








5 



0.50 

0.45 

0.40 

0.35 

0.30 

0.25 

0.20 

0.15 



0.50 

0.45 

0.40 

0.35 

0.30 

0.25 

0.20 

0.15 



0.40 

0.35 

0.30 

0.25 

0.20 

0.15 

0.10 

0.05 

0.00 

-0.05 



0.25 

0.20 

0.15 

0.10 

0.05 

0.00 


Figure 4. Short distance correlation functions {b^bx) (top panels) 
and {(no — \){nx — |)) (bottom panels) for the soft-core HFBI 
(left panels) and the hard-core FIFBI (right panels) on the FF = 8 
cylinder. The magnitude of the correlation function at site x is pro¬ 
portional to the radius of the corresponding circle. 


Soft-Core 
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Figure 5. Bound on the correlation length of all operators and 
(7(1) symmetric operators respectively for the soft-core and hard¬ 
core states vs. cylinder circumference W. Fits of the form ^ = 
were used to extract the correlation lengths. These 
bounds can be confirmed to match the correlation lengths of {blby) 
and (rixny) in each case. 


while the correlation functions in an MPS of finite bond di¬ 
mension decay exponentially. The long-range correlations of 
an MPS are encoded in its transfer operator T, which for an 
MPS of bond dimension M is a matrix of size x M^. De¬ 
noting the spectrum of T as Xi with |Ao| > |Ai| > ..., we 
can normalize the state such that Aq = 1. If the largest eigen¬ 
value is found to be non-degenerate, Ai < Aq, we have that 
all correlation functions of operators Oi that are supported 
on a finite number of sites centered around a site i decay as 
(OiOj) — {Oi){Oj) ~ Crucially, the correla¬ 

tion length fo for any operator O is bounded from above by 
— 1/ log I All [51]. In the following, we thus evaluate the spec¬ 
trum of the transfer operator of our PEPS along cylinders of 
varying circumference W to establish an upper bound on the 
correlation length for each circumference ^(W) Note that the 
possibility of having power-law correlations in a PEPS can 
be reconciled with the above consideration if the correlation 
length ^{W) diverges as FP —?► oo; we will thus need to care¬ 
fully consider the scaling of ^ {W). 

Our results for the correlation bounds ^(FP) are shown in 
Fig. 5. Here, we show the upper bound for the case of soft¬ 
core and hard-core bosons, and in each case consider the spec¬ 
trum of the full transfer operator as well as = 0 sector, 
which encodes correlations of operators Oi that do not change 
the boson number (such as density-density correlations). We 
find that in each case, the largest eigenvalue of the transfer 
operator is non-degenerate. Furthermore, we find that the cor¬ 
relation length approaches a finite constant as we increase FP, 
as shown in Figure 5. 


IV. ENTANGLEMENT SPECTRUM 

The quasi-ID cylinder geometry is convenient for calculat¬ 
ing the entanglement spectrum for entanglement cuts trans¬ 
verse to the long direction of the cylinder. Here, the entan¬ 
glement spectrum £i is defined through the spectrum pi of the 


reduced density matrices p^/n obtained for a bipartition of 
the state, where we have £i = — log pi. The corresponding 
eigenvectors of the reduced density matrices are referred to as 
Schmidt states The Schmidt decomposition 

IV’) = XI 

i 

relates the Schmidt states to the original wavefunction and is 
useful to keep in mind when interpreting the entanglement 
spectrum. 

To extract the entanglement spectrum exactly, we use a 
method proposed in Ref. 49. In the setup given here, the ex¬ 
act representation of the HFBI state as a PEPS of fixed bond 
dimension implies an upper bound on the number of non-zero 
pf, for the cut shown in Fig. 3(e), this upper bound is x = 2^. 

Upon computation of the entanglement spectrum, we find 
that this bound is saturated, so that there are precisely 2^ 
contributing terms in Eq. 6. This fact can be simply under¬ 
stood without reference to the PEPS representation: Each of 
the w plaquettes on the cut can contribute its one boson ei¬ 
ther to the left or the right of the cut, and this is the complete 
source of the uncertainty of the state on one half of the cut 
when ignoring the state on the other half. 

We can form a set of 2^ vectors |cti, ..., ctw) that cor¬ 
respond to the choices for the auxiliary degrees of freedom 
of the PEPS across the cut, where ai G {0,1} is the num¬ 
ber of bosons contributed by the Tth hexagon the left of the 
cut. The PEPS defines a map from these boundary vectors to 
physical states in the bulk of the semi-infinite cylinder, which 
we denote as on the subspace of physical states 

spanned by Schmidt states with non-vanishing contribution to 
the reduced density matrix, this map is invertible and can be 
computed explicitly. 

Translation around the cylinder acts in the natural way on 
the states permuting the values of ai, ai -G 
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K K 

Figure 6 . Entanglement spectrum on the zig-zag edge of a cylinder 
of circumference W — 10 (left panel) and W = 9 (right panel), 
as function of transverse momentum Ky. Different colors indicate 
different (7(1) charge sectors. 


CTi+i. Although the boson number on the half-infinite cylinder 
is infinite, we can define for each basis element a (7(1) charge 
cotTesponding to the number of bosons to the left of the cut 
relative to a uniform background charge, 

QL = ^{(yi-^)- (7) 

i 

Only relative charges between states will be important for our 
conclusions, and the precise way in which background charge 
is accounted for does not matter. Each state is paired with 
a corresponding state on the right half of the cylinder with 
opposite charge. 

We can block-diagonalize the reduced density matrix for 
a half-infinite cylinder in both the (7(1) charge and trans¬ 
verse momentum quantum numbers, allowing us to perform 
more efficient calculations. In addition, we can assign quan¬ 
tum numbers to both the Schmidt states and the entangle¬ 
ment spectrum. This property is generically true for (7(1)- 
symmetric and translationally invariant PEPS on a cylinder, 
although not in general true for arbitrary symmetry groups. 
This point is elaborated on in Appendix A. 

The entanglement spectra for the HFBI on cylinders with 
even and odd width circumferences are shown in Eig. 6, plot¬ 
ted against the transverse momentum eigenvalue and colored 
by the (7(1) charge eigenvalue of the corresponding Schmidt 
states. All the numerical results in this section are obtained for 
the soft-core boson variant of the state. We find that the entan¬ 
glement spectrum looks like it has a gapless edge mode with 
linear dispersion near momentum zero. To further substantiate 
this, we compare the lowest entanglement energies for several 
cylinder widths and quantum number sectors in Fig. 7. The 
hnite-size scaling confirms in all cases that the entanglement 
gap closes as 1/W, as you would expect for a gapless mode 
with linear dispersion. 

The gapless edge is suggestive of the state having either 
topological or SPT order. While topological order was al¬ 
ready ruled out in Ref. [32], our PEPS representation gives 
us additional tools to substantiate this assertion. In particular 
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Figure 7. The lowest five states above the ground state in Figure 6 
show gapless IjW scaling behavior. In this plot, fits for the entan¬ 
glement energy versus 1 /FK (of the form a + ^ + -^) to extract 
the gap are consistent with a gap value of 0 . 



Figure 8 . The constant correction to the entanglement entropy, which 
measures the topological entanglement entropy 7 when the state is 
minimally entangled, is consistent with 0 . 


we can calculate the entanglement between the two parts of 
the cylinder as a function of circumference S'(FP) and check 
for a subleading term to the area law by performing a fit to 
S'(FP) = aW + Sq. In a topological phase and in one of 
the minimally entangled states (MES) [52], one would ex¬ 
pect the subleading term to correspond to the topological en¬ 
tanglement entropy, Sq = —7 [53-55]. In a non-minimally 
entangled state, one would instead measure other values of 
^0 > — 7 . However, since each MES exhibits long-range or¬ 
der of a specific Wilson loop operator, such a superposition 
of MES represents a superposition of different ordering pat¬ 
terns and can thus be detected via a degeneracy of the largest 
eigenvalue of the transfer matrix. Our results for the entropy 
are shown in Fig. 8. We find results that are consistent with 
S'o = 0, which together with the fact that we also find that 
the transfer matrix to be nondegenerate rules out topological 
order. 
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A. Conformal field theory description of the edge 

In addition to the gapless behavior, we find that the low 
energy entanglement spectrum can be completely described 
by the finite-size spectrum of a conformal field theory (CFT). 
Given the C/(l) symmetry of the state, the simplest possible 
conformal field theory we might expect to appear at the edge is 
that of a single free bosonic field - and indeed, this is the CFT 
that matches the entanglement spectrum. We briefly review 
the relevant properties of this CFT [56]. 

The free boson CFT is created from the Lagrangian 

w 

£ = I J dt j dx (8) 

0 

with the compactified field identification 
(f) = (j) + 2 tiR 

and placed on a circle of circumference W with periodic 
boundary conditions 


(/>(x) = (j){x + W). 


The family of free-boson CFTs is parametrized by a single 
parameter k = irgR^, also known as the Luttinger liquid pa¬ 
rameter [56, 57]. 

Upon canonical quantization, we find that the set of energy 
eigenstates consists of t/(l) Kac-Moody primaries \e,m), 
with integers e,m labeling the [/(I) charge and the wind¬ 
ing number of the bosonic field respectively, and level n, n 
descendants of each primary for non-negative integers n, n, 
which we will collectively label |e, m; n, n). The number of 
level (n, n) descendants of a given primary, all of which are 
degenerate in the thermodynamic limit, is Z{n)Z{n), where 
Z[n) is the number of partitions of the integer n. 

The energies and momenta for the states \e,m;n,n) are 
given below on a finite size system of circumference W: 


^ 27r _ 

P = Trrlem + n — n] 

W 

27r 2 

H = 777 -1- Km +n + n) + ... 

W Ak 


(9) 


Here, the ellipsis (...) denotes further subleading contribu¬ 
tions due to coupling to irrelevant operators. 

By rescaling the energy and momentum, we find a system- 
size independent pattern that can be matched to the low- 
energy, linearly dispersing part of the entanglement spectrum 
from Figure 6: 


P oc {em + n — n) 

H cx -I- -I- 4K(n -|- n) -|- ... 


( 10 ) 


The results of this match are shown in Figure 9. An estimate 
for K can be obtained from the energy of the first descendent 
|0,0; 1, 0), which gives k « 1.6. The label e, which measures 
the U{1) charge, is integer for even W and half-integer for 



Figure 9. The identification of the primary states | ± e, m = 0) 
and the level n, n descendants in the spectrum of the soft-core boson 
entanglement Hamiltonian. The states are labeled e„,n. The zero 
and scale of the numerical spectrum are set by matching the lowest 
two states. The energies and charges of the primaries with charges 
2, 5/2, ...4 appear at the predicted spots. The best estimate for the 
Luttinger parameter from this spectrum is k « 1.6, taken from the 
energy of the Oi,o state. 


odd W. The degeneracy pattern 1,1, 2,... for the level- (n,0) 
descendents along the edge of the cone matches the predic¬ 
tion. 

The states with odd winding number m, such as |0,±1), 
do not appear at the energy and momentum predicted by the 
above formula. Instead, the primary states |e,m = ±1) can 
be found centered around momentum K = tt. The identifi¬ 
cation of these states in the spectrum is shown in Figure 10. 
Although the larger-m states are too high in energy to be re¬ 
liably distinguished at this system size, a natural conjecture is 
that all primaries with odd m will appear around momentum 
TT. (This is a standard side-effect of lattice regularization.) 

Given the PEPS representation, we can express the entan¬ 
glement Hamiltonian for the left semi-infinite cylinder, 
defined via pL — exp(—77^), as a Hamiltonian acting on 
the auxiliary degrees of freedom of the PEPS crossing the 
cut, which we have denoted as \ai,..., aw )• We expect this 
Hamiltonian to encode the universal properties of the edge 
GET, which should be invariant under local gauge choices in 
the PEPS. Its ground state is (up to normalization) given as 


\^o)= E (11) 

(7i,...,aw 


In Pig. 11, we show the bipartite von Neumann entanglement 
entropy of this state for a cut into I and W — I sites, which 
confirms the central charge c = 1 of the edge CPT A similar 
construction was considered in [58]. 
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Figure 10. The identification of primary states |e, m = ±1) and first 
descendants in the low energy part of the spectrum near momentum 
TT. Unlike the m = 0 states shown in Figure 9, these primary states 
have shifted momentum 7T = tt + em(2Ti/W) and an extra double 
degeneracy due to the two values of m = ±1. Using the estimate 
K « 1.6 from Figure 9, the predicted value of the entanglement en¬ 
ergy for the |e = 0, m = ±1) states is which has been marked 
in the plot. The agreement is very good. 



Number of sites in cut 


Figure 11. Entanglement entropy within the entanglement ground 
state of the soft-core boson state on 10 sites. For comparison, the 
Calabrese-Cardy formula [59] S{x) = c/31ogsin(7ra:/L) -f const. 
is shown with c = |, 1, and 2, with the const, fixed by matching the 
maximum of the entanglement entropy data, c = 1 is a good fit. 


V. SYMMETRY PROTECTION 
A. Overview 

While the gapless entanglement spectrum observed above 
is consistent with a symmetry-protected topological phase, it 
does not by itself guarantee the presence of such a robust 
phase, and does not allow us to infer which symmetries are 
protecting the topological properties of the phase. A key ob¬ 
servation that allows us to make progress on these crucial 
questions is that many points in the entanglement spectrum 
are degenerate. In particular, we find that for cylinders of odd 


circumference, the entire spectrum is doubly degenerate. In 
this section, we will discuss how the corresponding degen¬ 
erate Schmidt states are related through the action of a sym¬ 
metry of the HFBI wavefunction. As discussed in Ref. 45 and 
reviewed in the Appendix A, this symmetry action can be used 
to diagnose one-dimensional symmetry protected topological 
order, for which the degeneracy throughout the entire entan¬ 
glement spectrum is a robust feature. We will demonstrate 
that the odd circumference cylinders, considered as quasi-one- 
dimensional states, are indeed SPTs protected by a combina¬ 
tion of lattice inversion and charge parity symmetries. 

While the Schmidt eigenstates are uniquely defined for non¬ 
degenerate eigenvalues of the reduced density matrix, they are 
not unique when the spectrum is degenerate and any choice 
of orthonormal states in the degenerate subspace represents a 
valid choice of Schmidt states. Applying a unitary transfor¬ 
mation V^'', which respects = Sjk, on the left 

Schmidt states must be accompanied by an appropriate trans¬ 
formation (U-t*)* applied to the right Schmidt states. 

In particular, this allows the action of an on-site symmetry 
(or more generally, any symmetry which commutes separately 
with the reduced density matrices for the left and right half) to 
mix Schmidt states corresponding to degenerate eigenvalues. 
The action of such a symmetry operator Ug takes the form 

Ug\^L) = E 

^ ( 12 ) 

k 

where the are unitary matrices that only act on degenerate 
blocks of Schmidt states. Crucially, Ref. 45 describes a nu¬ 
merical procedure to calculate Vg for an on-site symmetry g 
within the MRS formalism, which we review in Appendix A. 

We can also analyze the effects of symmetries that preserve 
the entanglement cut but swap the left and right halves of the 
cylinders. In general, we will consider any symmetry h that 
swaps the cylinder sides and squares to the identity, which we 
will call an inverting symmetry. These satisfy a modification 
ofEq. (12): 


Uh\^f) = E 

^ (13) 

k 

Note here that the left and right Schmidt states are exchanged 
in the transformation. We can introduce a map S that acts as 

= \4^)- (14) 

Since a change in phase —>■ must be accom¬ 
panied by the complex conjugate —>■ to pre¬ 

serve the Schmidt decomposition, S is antiunitary. 

Combining the above, we see that 

UhSlPj^^) = E (15) 

3 
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defines the action of the operator UhS on the right Schmidt 
states (of course an equivalent action can be defined on the left 
Schmidt states). Since S is anti-unitary, the combined action 
of UhS is also anti-unitary. Together with the requirement that 
the symmetry squares to the identity, one finds that (where K 
represents complex conjugation in the canonical basis) 

VhV: = {VhKf = = ±1, (16) 

that is the inverting symmetry forms an anti-unitary projective 
representation of Z 2 . 

As reviewed in Appendix A, the collection of Vg for on¬ 
site symmetries sometimes fail to satisfy the group multipli¬ 
cation laws, i.e. one may find 7 ^ ^ 1 ^ 92 - Instead, they 
may form a projective representation, where group multipli¬ 
cation laws are obeyed up to phases uj{gi,g 2 ), i.e. Vg^Vg^ = 
92)Vgig2- Certain combinations of these phases, such as 

e* 0«„.2 ^ ^[ 91 , 92 ) 

^( 92 , 91 ) 

whenever [gi , ( 72 ] = Oj symmetry protected topological 
invariants, which take discrete values and hence cannot be 
changed continuously. Thus, 7 ^ 0 signifies that the 

entanglement degeneracy cannot be removed without break¬ 
ing the symmetry or going through a phase transition. Simi¬ 
larly, for the inverting (anti-unitary) symmetries h, the phase 
(ph = TT in Eq. (16) signifies that the degeneracy cannot be 
removed without breaking the symmetry [45]. 


B. Symmetry protection of the HFBI 

The on-site symmetries of the featureless boson insulator 
considered here are the U{1) charge symmetry and the anti- 
unitary time-reversal symmetry r, which acts by complex 
conjugation in the boson number basis. Despite being at half¬ 
filling, the hard-core boson variant of the state does not have 
a particle-hole symmetry. Exploring the edge action of these 
symmetries numerically, we find that they are all represented 
linearly and thus do not protect the degeneracy of the entan¬ 
glement spectrum on cylinders of odd circumference. In order 
to protect the degeneracy, we must therefore include lattice 
symmetries. 

By choosing a cylinder geometry, we explicitly break some 
of the lattice rotational and reflection symmetries. The re¬ 
maining symmetries are generated by translations parallel 
and Ty perpendicular to the cylinder axis as well as reflec¬ 
tions Ir about a line parallel and ly about a line perpendic¬ 
ular to the cylinder axis. We also consider lattice inversion 
I = Ixly, equivalent to a tt rotation of the spatial plane about 
the center of a hexagonal plaquette. These symmetries are 
illustrated in Eig. 12. We find that a number of symmetry- 
protected topological invariants that are defined through these 
symmetries take non-trivial values in the HEBI, thus protect¬ 
ing the doubly degenerate entanglement spectrum on odd cir¬ 
cumference cylinders. The complete list of non-trivial invari¬ 
ants is summarized in Table I. 


/ ^ ^ y 

\ _ / \ 
/ \ _ / 
\ _ / \ 
/ V^\ X, / 


/ v 
\ _ / 
/ \ 
\ _ / 
/ \ 


C3 




/ 

' \ 


Figure 12. Lattice symmetries considered here: (i) reflection 
about a line parallel to the long direction of the cylinder, (ii) Ty re¬ 
flection about a line perpendicular to the long direction, correspond¬ 
ing to the entanglement cut shown in Fig. 3(e). These are both chosen 
such that the reflection line crosses the hexagon center. Their prod¬ 
uct, T = TTTy, thus represents (iii) the inversion about a hexagon 
center. 


The crucial ingredient underlying these SET invariants is 
a spatial symmetry h that swaps the two sides of the entan¬ 
glement cut. By a general symmetry analysis of Eq. (15), 
Vh must act as a particle-hole symmetry on the edge, since 
the Schmidt pairing S (Eq. (14)) pairs states with opposite 
quantum numbers. In this case, the symmetry action Vxy is 
precisely that of a particle-hole transformation in the local 
PEPS basis. Defining \a) — |cri,..., crw) and |1 — f?) = 
|1 — (Ti, ..., 1 — (Tw), we And that 

Fija) = 11 -ct), (18) 

since a state where the hexagon contributes ai bosons on 
the right is paired with a state where the hexagon con¬ 
tributes 1 — (Ti on the left. We can thus read off that Vi^ acts 
like (Jx in the space spanned by the states {|(t) , |1 — (?)}. 

When W is odd, these states have opposite charge parity. 
Specifically, if II = G C/(l) is the charge parity sym¬ 
metry, we have 

Vz7|(t) = (-1)^-M?) 

= (_1)VF(_i)E(7.|i-^) (19) 

Therefore, for W odd, V 77 acts like cr^ in the space 
{jf?), |1 — (?)}. It is thus reasonable to expect that the combi¬ 
nation of these two symmetries acts as Vnx = UxCTz, which 
would obey the property that VnxVnx = *^I**^® form a 

topological invariant. 

The local PEPS basis is not unitarily equivalent to the 
canonical form basis, so we must check this numerically by 
performing an explicit calculation in the canonical form of an 
MPS representation of the state, as outlined in Appendix A. 
We thus confirm SPT invariants for symmetries that involve 
such a spatial symmetry h and an on-site symmetry. There are 
several appropriate invariants, as listed in Table I; the simplest 
is 


VnxV^x = -I- ( 20 ) 
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Group 

Generators 

Invariant 

i 


{nx} 

VnxV^i = 

-I 

— 


{nxy} 

VniyV^x, = 

-I 

- 

Z 2 X Zf ^ 

{77, rZ} 

VnVrxV-^V-^ = 

-I 

+ 

Z 2 X Zl’^ 

{n.rXy} 

VnVrXyV^^V-^\ = 

-I 

+ 

Z 2 X Zf 

{nx^,TX} 

Vnx^V^xV-lV-^^ = 

-I 

+ 

Z 2 X Zf 

{IIXx,TXy} 

Vnx^VrXyV-lV-^\ = 

-I 

+ 


Table I. Summary of symmetry protecting invariants found for the 
HFBI state. The last column indicates whether the symmetry acts 
unitarily {i = +) or antiunitarily {i = —1) on the edge. The de¬ 
generate entanglement spectrum cannot be split unless all 6 of the 
minimal protecting symmetry groups are broken. 


From this we see that the charge, translation, and inversion 
symmetry can all be broken without splitting the entanglement 
degeneracy, as long as the single combined symmetry IIX is 
preserved. In Section VI, we will discuss perturbations that 
preserve this symmetry. 

We note that there are also symmetries that act unitarily on 
the edge and yield SPT invariants; however, these must form 
the group Z 2 x Z 2 as Z 2 does not have unitary projective rep¬ 
resentations. Examples for this are formed by involving time- 
reversal symmetry; since I 4 and Vx both act antiunitarily, Vri 
acts unitarily on the edge. The Z 2 x Z 2 group generated by 
tX and 7T has a projective representation characterized by the 
topological invariant 

VnVriVn^V-^^ = -I. (21) 

This symmetry protection gives a distinct class of perturba¬ 
tions that cannot split the entanglement degeneracy. The com¬ 
plete set of symmetry groups we find is summarized in Table I. 

We can form variants of the HFBI state, which are unitarily 
related to the original state by an on-site unitary and thus share 
the same entanglement spectrum, where the entanglement de¬ 
generacy can be protected by a lattice symmetry alone without 
involving the on-site 77 symmetry. Essentially, the protecting 
symmetries of the variant generated by a unitary U is obtained 
by conjugating the generators of the protecting symmetries of 
the HFBI by U. These will be discused further in Appendix C. 

We also mention that the symmetry protected invariants 
produced above imply the existance of non-local correlations 
in the form of ‘membrane’ order parameters that naturally 
generalize the string order parameters from one-dimensional 
SPT phases [60]. For example, the sign of V /71 V^x 
tected by measuring the overlap of the state |'0) and the same 
state with a partial application of the protecting symmetry, i.e. 

lim {'ijj\{nx)i 2 n\i’) OC (-1)^, (22) 

n—)-oo 

where (TTIji ^n is the restriction of UX to 2n cylinder slices. 
We leave open the question of whether this ‘membrane’ order 
parameter generalizes in any way to regions that do not wrap 
the cylinder. 


C. Tight-Binding Restriction 

Before concluding this section, we comment briefly on the 
role played by the restriction to a particular tight-binding 
model. Restricting to a particular tight-binding model is a 
stronger condition than merely specifying a space group sym¬ 
metry. For instance, the triangular, honeycomb, and kagome 
lattices all share the same space group, but encode it using 
one, two and three orbitals per unit cell, respectively. An ex¬ 
ample is graphene: the electronic Dirac cones in graphene 
are protected (for vanishingly small spin-orbit coupling) not 
solely by its space group symmetries, but rather by the tight- 
binding representation of those symmetries [61] [62]. A re¬ 
striction to a tight-binding representation is often well mo¬ 
tivated by experiments. Choosing a tight-binding model 
amounts to restricting to the class of models that can be rep¬ 
resented using precisely the orbitals we began with; if we are 
given the freedom to add sites, it may be possible to exit a 
topological phase and enter a trivial one without closing the 
gap, while still preserving lattice symmetry. For the HFBI, 
this would be accomplished by adding sites at hexagon cen¬ 
ters, and adiabatically deforming the relative weights afforded 
to bosons placed at the new sites and the original honeycomb 
lattice positions. The symmetry protection of the entangle¬ 
ment structure is thus fairly subtle. 

This subtlety is better understood for the simpler case of 
non-interacting fermions. The classification of free-fermion 
topological phases is known to be richer if one removes the 
freedom to add trivial bands [63]. The continuity between two 
phases which arises upon the addition of such trivial bands 
is known as “stable equivalence” in accord with a basic no¬ 
tion in K-theory [63]. We note that this tight-binding re¬ 
striction also distinguishes the HFBI and related symmetry- 
protected short-range entangled states from the ground states 
of “filling-enforced” topological band insulators introduced 
very recently [64] . 


VI. QUASI-LOCAL PARENT HAMILTONIAN AND 
PERTURBATIONS 

We now re-examine the question of whether the HFBI state 
is representative of a robust phase of matter that is separated 
from conventional phases by phase transitions. One way to 
demonstrate this - beyond the topological invariants discussed 
above - is to find a local Hamiltonian with a unique ground 
state that is the HFBI wavefunction and study the ground state 
properties under perturbations to this Hamiltonian. For many 
tensor network states - those that satisfy an injectivity condi¬ 
tion [65] — a frustration-free, local parent Hamiltonian with a 
unique ground state can be explicitly constructed. In our case, 
this injectivity condition can be shown to not hold on any sim¬ 
ply connected cluster of sites that we can numerically access, 
and this specific construction of a parent Hamiltonian is thus 
not possible. Given the challenges of numerical simulations 
of two-dimensional systems, an exhaustive numerical search 
for such a Hamiltonian seems unfeasible. 

To avoid these problems, we will focus on a quasi-ID ap- 
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proach in this section. This is based on the observation that 
while the PEPS is not injective on simply connected clusters, 
it does turn out to be injective on slices of an infinite cylinder. 
This gives rise to a gapped Hamiltonian whose unique ground 
state is the HFBI. This ‘parent Hamiltonian’ is local in the 
non-compact direction of the cylinder, but non-local around 
the cylinder and dependent on the circumference W. We be¬ 
lieve that nevertheless, the insights gained from these (par¬ 
tially non-local) Hamiltonians can serve as a starting point for 
identifying the phase in more sophisticated numerical studies 
of fully two-dimensional boson systems. 

Given the unperturbed Hamiltonian, we study the robust¬ 
ness of the entanglement spectrum to perturbations. This de¬ 
pends on the class of perturbations allowed - SPT phases are 
only distinct if perturbations that break the symmetry are for¬ 
bidden, which is reflected in the fact that the topological in¬ 
variants that distinguish the phases are ill-defined in the ab¬ 
sence of the symmetry. According to the results discussed in 
Section V, it is not necessary for the perturbations to preserve 
the entire symmetry group of the HFBI wavefunction to pre¬ 
serve the entanglement in the state. Instead, the entanglement 
is robust to any perturbation that does not break all of the six 
protecting groups discussed in Table I. This set of perturba¬ 
tions is much bigger than the set of perturbations that preserve 
the entire symmetry group of the HFBI wavefunction. We will 
confirm that the double degeneracy throughout the entire en¬ 
tanglement spectrum survives these perturbations for odd-W 
cylinders, while it splits for other perturbations that break all 
of the protecting symmetries. 


A. Parent Hamiltonians for the W = 1 cylinder and 
equivalence to the Haldane insulator 

The W = 1 cylinder with hard-core bosons has a Hilbert 
space equivalent to a two-leg spin-^ ladder. The HFBI state in 
this case has a natural MBS representation of bond dimension 
d = 2, constructed by contracting the tensors around each 
cylinder slice. A well-known property of MBS is the existence 
of a parent Hamiltonian - a frustration-free Hamiltonian with 
its unique ground state given by the MBS, first introduced by 
Ref. 38. The parent Hamiltonian is constructed in this case 
as a translationally invariant sum of projectors, where each 
term projects the Hilbert space of two consecutive rungs of 
the ladder to the = 4 dimensional subspace of states form 
spanned by the nonzero eigenvectors of the reduced density 
matrix of those two rungs. The result Hq of this construction 
involves all possible terms that act on two rungs of the ladder 
and preserve charge and reflection symmetry. 

Using a local unitary transformation discussed in detail in 
Appendix B, we can transform the wavefunction of the HFBI 
on the W = 1 cylinder to that of the ‘Haldane insulator’ [21, 
45], which is known to be the ground state of an extended 
Bose-Hubbard model on the two-leg ladder in an appropriate 
parameter regime, and has also been shown to be a ID SBT 
with a doubly degenerate entanglement spectrum and a non¬ 
local string order parameter protected by charge parity times 
inversion III. 


This extended Bose-Hubbard Hamiltonian that gives rise 
to the Haldane insulator includes only hopping and density- 
density interactions; additionally, the range of these interac¬ 
tions extends only to neighboring rungs of the ladder. It is 
thus clear that the additional interactions present in the parent 
Hamiltonian for the HFBI can be tuned away without under¬ 
going a phase transition. The hard-core bosons should addi¬ 
tionally be considered to have infinite on-site density interac¬ 
tions, which can be tuned away from infinity to make a state 
with soft-core bosons. 


B. Perturbing the state on the IT = 3 cylinder 

Similar to the W = 1 cylinder, we can obtain a parent 
Hamiltonian for the W = 3 cylinder as a sum of local projec¬ 
tors acting on adjacent slices of the cylinder. We then consider 
two different perturbations to these quasi-local parent Hamil¬ 
tonians. For each perturbation, we use infinite time-evolving 
block decimation (iTEBD) [66-69] to evolve an initial wave- 
function in imaginary time until it converges to the ground 
state of the perturbed Hamiltonian. The two perturbations 
considered are the superfluid pairing 

H' = AY,hb,+h.c., (23) 

(U> 

which breaks the (7(1) charge symmetry down to the Z 2 
charge parity subgroup and the uniform held 

H' = hY,{h + bl), (24) 

i 

which fully breaks (7(1) charge symmetry but preserves lat¬ 
tice symmetry. The perturbation in (23) does not break 
the protecting symmetry III, while the perturbation in (24) 
breaks all of the protecting symmetries. 

Figure 13 show the resulting entanglement spectra from the 
ground states obtained with iTEBD. Indeed, the perturbation 
in (24) splits the degenerate entanglement spectrum, whereas 
the double-degeneracy of the entire spectrum is preserved for 
those perturbations that do not break all of the protecting sym¬ 
metries. In the case of a symmetry-breaking perturbation, the 
splitting is most easily observed for the higher levels, but - 
as shown in the inset - also the lowest values of the entan¬ 
glement spectrum are weakly split by an amount that scales 
roughly linearly in the strength of the perturbation. 

Unfortunately it is beyond the scope of this work to de¬ 
termine which perturbations leave the CFT structure of the 
entanglement spectra intact. To assess this would require 
a Hamiltonian that is local in two dimensions (rather than 
the Hamiltonians used here which are only local in the non¬ 
compact direction of the cylinder). Furthermore, it would 
require being able to perform accurate simulations for large 
cylinders, which is prohibitive with the techniques used here. 
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Figure 13. Entanglement spectrum for the ground state of the par¬ 
ent Hamiltonian on the W = 3 cylinder under a symmetry-breaking 
perturbation (24) (left panel) and a symmetry-preserving perturba¬ 
tion (23) (right panel). These results obtained using an iTEBD sim¬ 
ulation using a bond dimension of M — 24. 


VII. BOSON-VORTEX DUALITY 

An alternative approach for examining the phase structure 
of two dimensional bosonic systems is to use the boson-vortex 
duality, which rewrites the theory in terms of the superfluid 
vortex degrees of freedom defined on the dual lattice. In 
this picture, the site filling of bosons on the original lattice 
is mapped to an effective magnetic flux through dual lat¬ 
tice plaquettes that modifies the vortex hopping via the usual 
Aharonov-Bohm phases [70, 71]. In the dual description the 
superfluid and Mott insulating phases of the bosons are re¬ 
spectively mapped into the gapped and condensed phases of 
the vortices. It is instructive to see how this approach fares on 
the honeycomb lattice at half-site-filling. The vortices move 
on the dual triangular lattice, and the original site filling of 
1 /2 corresponds to a tt-Aux for vortices for every triangular 
lattice plaquette. Each unit cell on the triangular lattice con¬ 
tains a pair of triangles and hence 27r flux; as a consequence, 
the vortices transform normally (rather than projectively) un¬ 
der lattice symmetries. Naively, the tt-Hux has the effect of 
inverting the vortex band structure so that the vortex minima 
are shifted to the Brillouin zone corners K, K' rather than the 
zone center T. Condensing vortices at the zone corners would 
break lattice symmetries [72]. However, the fact that the vor¬ 
tices transform regularly under lattice symmetries — in other 
words, that the flux pattern does not lead to an enlarged mag¬ 
netic unit cell — allows us to add additional hopping while 
preserving symmetries to returns the vortex minimum to E. 
Condensing vortices at F then restores the (7(1) symmetry 
while preserving all lattice symmetries. As there are no other 
known symmetry-preserving insulating phases of bosons on 
the honeycomb lattice at this Ailing, we conjecture that this 
phase is adiabatically connected to the HFBI. 

It is perhaps worth noting that the vortex-condensation pic¬ 
ture also illustrates a fundamental distinction between half- 
filling on Bravais and symmorphic non-Bravais lattices. As an 
example, consider half-filling on the square lattice [73]; per¬ 
forming the duality transformation, we arrive at a theory of 
vortices moving on the dual square lattice with tt flux through 
each plaquette. Crucially, this flux assignment on the square 


lattice doubles the unit cell, and so vortices form a projec¬ 
tive representation of the space group (related to the mag¬ 
netic translations familiar from studying particles in a mag¬ 
netic held). This projective structure cannot be removed and 
so guarantees that a single non-degenerate minimum cannot 
be restored for any choice of vortex hopping parameters. Put 
differently, the vortices are (unlike in the honeycomb case) 
forced to carry non-trivial space group quantum numbers, 
and the condensation of single vortices necessarily breaks the 
symmetry [73]. Other approaches lead to more exotic alter¬ 
natives, e.g. condensing vortices in pairs triggers fractional- 
ization. This in accord with the expectation that a gapped fea¬ 
tureless insulator is absent on the half-filled square lattice [26- 
29]. 


VIII. CONCLUDING REMARKS AND DISCUSSION 

We have applied recently developed tensor network meth¬ 
ods to study the edge properties of a bosonic insulator that is 
featureless in the bulk. Our simulations are performed for an 
infinitely long cylinder of finite circumference W. This allows 
us to numerically extract the exact entanglement spectrum for 
up to FF = 10. We And that the entanglement gap closes as 
1/FF, and that furthermore the low-lying spectrum coincides 
to high accuracy with the spectrum of a free boson confor¬ 
mal held theory. This is further corroborated by observing a 
central charge of c = 1 in the entropy of the lowest Schmidt 
state. 

While these observations are consistent with and strongly 
suggestive of a symmetry-protected topological phase, where 
such a gapless spectrum would naturally emerge at the edge, 
these calculations do not establish a rigorous connection be¬ 
tween the edge spectrum and symmetry-protection, i.e. they 
leave open the possibility that the gapless entanglement spec¬ 
trum is accidental. To make progress on this question, we ana¬ 
lyze in some detail the exact degeneracies in the entanglement 
spectrum for cylinders of odd circumference W. Using re¬ 
cently developed tools based on matrix-product states, we are 
able to establish a strong connection to the symmetries of the 
state by computing topological invariants that detect the non¬ 
trivial action of certain symmetry operations. These symmetry 
operations, whose action is non-trivial, consist of particular 
combinations of lattice and spin symmetries on the edge. This 
establishes in the affirmative that the quasi-one-dimensional 
systems obtained for odd cylinder widths W represent one¬ 
dimensional symmetry-protected topological phases. 

We cannot establish with the same rigor that the symme¬ 
tries that protect these one-dimensional topological invariants 
also protect the gapless edge spectrum on the edge of the two- 
dimensional system. However, several considerations are in 
favor of this. Firstly, we observe that the symmetries that are 
shown to be relevant to the case of odd W are not inherently 
one-dimensional and could apply equally well to the full, two- 
dimensional system. The partial application of symmetry in 
the non-local order parameter in Eq. (22) could be applied 
to arbitrary inversion-symmetric regions in the plane, and not 
only to cylinder slices. Additionally, we can construct an ar- 
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gument based on the picture of the edge physics provided by 
the tensor network representation. As outlined in Section IV, 
the edge of the tensor network representation with the cut cho¬ 
sen here can be represented using the Hilbert space of a model 
of hard-core bosons hopping on a one-dimensional chain with 
one site per plaquette, where the occupation of a site corre¬ 
sponds to whether the boson of that plaquette is found on the 
left or right side of the cut. In this representation, the reflec¬ 
tion symmetry about the cut takes the special role of guar¬ 
anteeing equal probability for the boson to be on the left or 
right of the cut, and thus fixing the model for the edge to half- 
filling. Thus, if the edge remains translationally symmetric, 
our model for the edge has fractional charge per unit cell. If 
the entanglement Hamiltonian can be thought of as local, the 
Lieb-Schultz-Matthis theorem applies and guarantees that the 
entanglement edge is either gapless or spontaneously breaks a 
symmetry. This suggests that the phase is a two-dimensional 
symmetry-protected topological phase with a protecting group 
that includes translation and UX. 

The calculations presented here provide a case study where 
tensor network representations lead to novel insights into 
strongly correlated physics beyond what is accessible to more 
traditional methods, such as the quantum-to-classical map¬ 
pings pursued in Ref. 32 and reviewed in Section II. The 
tensor-network techniques used in the present approach allow 
us to strengthen the conclusions of Ref. 32, in particular on the 
absence of topological order, and reveal entanglement proper¬ 
ties that are entirely out of reach of quantum-to-classical map¬ 
pings. It is amusing to note that this development in theoreti¬ 
cal methodology closely parallels the history of the prototypi¬ 
cal SPT phase, the AKLT phase of the spin-1 chain, where the 
existence of a quantum-to-classical mapping was known well 
before the nontrivial entanglement structure was understood. 
As in that example, we expect that here as well, the quantum- 
to-classical mappings are restricted to rather special points 
within a broader SPT phase, whereas the tensor-network de¬ 
scription and its corresponding entanglement structure are ex¬ 
pected to be valid more generally throughout the phase. 

The question of a parent Hamiltonian, i.e. whether the 
HFBI can be established as the unique ground state of a 
gapped local Hamiltonian, remains open. As reviewed briefly 
in Sec. VI, the structure of the PEPS does not allow us to 
straightforwardly extract a local parent Hamiltonian in two di¬ 
mensions. However, this by no means implies that such a par¬ 
ent Hamiltonian does not exist, and future work will explore 
different numerical approaches to And such a Hamiltonian. 
Note: While completing this work, we became aware of re¬ 
lated PEPS constructions of featureless paramagnetic wave- 
functions on the square lattice with spin 1 per site, and on the 
honeycomb lattice with spin 1 or ^ per site [74]. The spin- 
i honeycomb lattice example corresponds to the same Ailing 
as the featureless insulating phase considered here, but has 
higher symmetry {SO{3)) compared to the 17(1) symmetry in 
the present paper. We note that even in the case where we con¬ 
sider spinful fermions (see Appendix C 2) bound into Cooper 
pairs, the wavefunction we construct here is not a valid wave- 
function for a spin-only model: projecting it to the case of 
single-fermion occupancy per site (as appropriate to a spin 


model) annihilates the wavefunction. It will be interesting to 
study if the spin-only wavefunctions constructed in Ref 74 
possess similarly rich entanglement structure as the HEBI. 
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Appendix A: Determining the edge action of the symmetry 
using MPS 

We can use the formalism of matrix-product states to deter¬ 
mine the action of physical symmetries on the Schmidt states. 
Eirst, this will lead to the assignment of charge and translation 
(which both act on-site in the MPS representation) quantum 
numbers to the Schmidt states and corresponding entangle¬ 
ment spectrum as labeled in, e.g.. Pig. 6. Secondly, this will 
be used to numerically extract the topological invariants dis¬ 
cussed in Section V. We now review this formalism briefly, 
including a discussion of the method that allows us to nu¬ 
merically determine the symmetry action of inversion sym¬ 
metry on the Schmidt states. Both of these discussions follow 
Ref 45. 

We start by finding tensors T, A representing the so-called 
canonical form of the MPS, as detailed in Refs. 66 and 68: 

\tjj) = ^. ..ATp^ATp^ATp^A... |...poPiP2---)- (Al) 
iPi} 

This canonical form provides the Schmidt decomposition at 
each site in the lattice. Here, each physical leg of the MPS 
represents all 2W physical sites on a cylinder slice, and each 
virtual leg represents all virtual indices that connect cylinder 
slices; the bond dimension of the MPS is thus 2'^. The change 
of basis to canonical form generally mixes the Hilbert spaces 
from these virtual legs, so the resulting basis won’t be local 
around the circumference of the cylinder. 

Por each on-site symmetry of the wavefunction Ug = 
with C/glV') = \'tp), there is an operator Vg that acts on the 

virtual leg of the MPS and satisfies the equation 





This equation can be rewritten and solved as an eigenvector 
problem; for an MPS with a nondegenerate largest transfer 
matrix eigenvalue, this equation is guaranteed have a unique 
solution where the eigenvalue is the largest eigenvalue 
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of the eigenvector problem. These solutions Vg have two im¬ 
portant properties: they are only defined up to a phase, and 
they are guaranteed to commute with the diagonal matrix A of 
Schmidt weights. 

Due to the first property, these operators are not guaranteed 
to obey the group multiplication laws, i.e. one could find situ¬ 
ations where 

VgVh=u:{g,h)Vgh. (A3) 

It is not always possible to absorb these phases into the deh- 
nitions of the Vg; in those cases, the Vg do not form a linear 
representation of the group but rather a projective represen¬ 
tation. The set of equivalent classes of phases uj{g, h) under 
redefinitions Vg —>■ a{g)Vg is called H^{G, U{1)), the second 
group cohomology with C/(l) coefficients. 

For all the groups discussed in this paper, the group coho¬ 
mology classes are labeled by elements of a discrete abelian 
group. These discrete classes cannot be connected to each 
other continuously without undergoing a bulk phase transi¬ 
tion or breaking the symmetry. Additionally, the classification 
of projective representations for the on-site symmetry group 
?7(1) X Z 14 / representing charge and translation around the 
cylinder is trivial. Thus, these edge symmetries can be taken 
to act linearly, and all Schmidt states can always be simul¬ 
taneously assigned charge and momentum eigenvalues, as in 
Figure 6 . 

The second property guarantees that the Vg only mixes ex¬ 
actly degenerate Schmidt states. The action of Vg must have 
the same phases uj{g,h) on each degenerate block of Schmidt 
states, so the projective representation can be nontrivial on 
any block only if every Schmidt state throughout the entire 
spectrum is degenerate. The degeneracy will be protected by 
the symmetry if and only if the Vg form a nontrivial projec¬ 
tive representation. Therefore this ID SPT analysis can only 
potentially give a nontrival answer for the odd W states of 
the HFBI, where this exact degeneracy is seen throughout the 
spectrum. 

The MPS analysis of inversion symmetry proceeds simi¬ 
larly. We will consider in general any symmetry h of the 
wavefunction that squares to the identity and that can be writ¬ 
ten in the MPS as the product of an on-site symmetry action 
Uh and a transpose of the site tensor. This will include an 
inversion of the honeycomb lattice - equivalent to a 180 de¬ 
gree rotation about the center of any plaquette, which we label 
I — lyXx, and the combination of inversion with on-site sym¬ 
metries. In addition, by blocking two site-tensors together, we 
can write the reflection symmetry Xy in this form as well. In 
this scenario, the edge symmetry action satisfies 



The map Vh is also computed from an eigenvector problem. 

For the HFBI, the symmetry group respected by the cylin¬ 
der geometry is 17(1) x x Z 2 ) x Z^ x Z|’, where the 


factors refer to charge symmetry, translation around the cylin¬ 
der, Xx, Xy, and r respectively. The P and T denote space- 
reversing and time-reversing symmetries, and signify the an¬ 
tiunitary action on the Schmidt states. Many of the non-trivial 
projective representations of such a complicated group will re¬ 
main projective when the symmetry is restricted to a subgroup 
- in this case, the full symmetry is not needed to protect the 
entanglement degeneracy. As shown in Table I, the projective 
representation corresponding to the HFBI state can indeed be 
protected by any one of a number of subgroups of the full 
symmetry group, all involving inversions and charge parity. 

The symmetry actions - both on-site and inversion symme¬ 
tries - are computed in the Schmidt basis, but can be trans¬ 
formed into the basis |{cri}) determined by the virtual legs of 
the PEPS in Figure 3(e). In this case, the symmetry action Vi^ 
is precisely a particle-hole symmetry in the local PEPS basis, 
with coefficients 

|o’i,...,crw) = \1 — ai,... ,1 — aw), 

since a state where the hexagon contributes ai bosons on 
the right is paired with a state where the hexagon con¬ 
tributes 1 — (Ti on the left. Thus 

i 

where K is complex conjugation in the local PEPS basis, and 
af is the Pauli operator acting on the site of the local PEPS 
basis. 

Charge symmetry acts locally as well: 

1 ^^,..., 1 ^^, ...,aw). 

In particular, charge parity Vn = e*®® can be written as 
Vn = 

i 

The combined action of charge parity and reflection across 
the cut takes the form 

V^Iy=\{{^a\)K, 

i 

which is precisely the form that time-reversal acting on an 
ordinary spin-7 chain takes. When the circumference of the 
cylinder W is odd, we see that 

VnXyVj^Xy = -I- 

The degeneracy of the entanglement spectrum can be seen as 
an application of Kramer’s theorem. Eormally, this property 
is said to characterize the nontrivial projective representation 

i72(Zf;[/(!)) =Z2, 

and remains true while XlXy is a symmetry and no phase tran¬ 
sitions have occurred. 





15 


Time reversal symmetry acts as complex conjugation in the 
local PEPS basis 14 = K. Translation and act as permuta¬ 
tions of the local PEPS basis: 

Vt\(Jit ■ ■ T<Jw) = |o'2,--.,crvi/,cri) 

Vl^ |(Ti,...,crw) = 


Each spin-1 degree of freedom can be split into a pair of 
5 = 1/2 spins to make a state on a spin-i ladder. An appro¬ 
priate Hamiltonian can be found in terms of the new spin vari¬ 
ables Sj^A/B by adding a term to project out the spin-singlet 
component of Sj^a + Sj^B- The spin-^’s can then be treated 
as hard-core bosons. The Hamiltonian becomes 


These symmetries can be combined with to create the 

additional topological invariants shown in Table I. A non¬ 
trivial projective representation in 

H \ Z 2 X Z2;C/(1)) =Z2 

is created whenever two unitary symmetries that commute in 
the bulk satisfy 


Va.v,,v-X-^ = -i. 

Each new invariant is related to a new set of pertubations that 
can’t break the entanglement degeneracy. 


Appendix B: From the AKLT to the W = 1 HFBI 

The AKLT state |'!/’aklt) is a state of a spin-1 chain that 
has an exact representation as an MPS of bond dimension 2 
using site tensors A? related to the Pauli matrices [40]. It is 
the exact ground state of the AKLT Hamiltonian 

TfAKLT = X/ + l + 

j 

but it is known that the simpler Hamiltonian 

Haf = Sj+i (B2) 

3 

is in the same phase, i.e. the AKLT state lies in the Haldane 
phase of the spin-1 Heisenberg chain. By a series of trans¬ 
formations, we can find a representative MPS wavefunction 
|4hi) and a simple representative Hamiltonian that can be adi- 
abatically connected to the W = 1 HEBI and its correspond¬ 
ing parent Hamiltonian. 

By using the unitary operator 

U{-k) = (B3) 

j even 

which flips the cc, y components of the spins on every other 
site, we create a wavefunction representative of the Haldane 
insulator (HI) [21] phase, which is protected by UlW — 
nX [45]. This phase is obtained as the ground state of the 
Hamiltonian 

H' = UHafU^ (B4) 

= E ■ (B 5 ) 

3 ^ ^ 


^ “2 + h.c.) 

3 

+ V{njA + rijB - l){nj+iA + nj+iB - 1 ) 

- + h.c.) - J{njA - - ^), (B6) 


where f = 1, L = 1, and J —> c». The J term projects the 
spin-singlet out of each rung, and in practice only needs to be 
larger than all other relevant scales to drive the system into the 
appropriate phase. 

We can do the same transformations on the MPS IV'aklt) 
to obtain a new MPS IV'hi) with bond dimension 2 and site 
tensor AT that represents a state in the phase of Hbi on the 
two-leg ladder. The site tensor Sf^ of the W = 1 HEBI 
also has bond dimension 2 and represents a state of hard-core 
bosons on the two-leg ladder. Numerically, these are repre¬ 
sented by the (unnormalized) site tensors 


and 



p = (00) 

P = (01) 

p=(10) 
p = (11) 

p = (00) 
p = (01) 
p = (10) 

p = (11), 


where p = {P 1 P 2 ) represents the occupation numbers of the 
hard-core bosons on the two sites on each leg of the ladder. 
By linearly tuning the site tensors using 

+ (B7) 

and checking that the transfer matrix of the resulting state is 
non-degenerate for all t G [0,1], we confirmed that the IL = 1 
HEBI can be tuned in the space of bond dimension d = 2 
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MPS to l^/’Hi) without passing through a phase transition, and 
the representative Hamiltonian in (B6) describes a state in the 
same phase. 

By calculating the canonical form of the Sf^ site tensor, one 
can check that the W — 1 HFBI wavefunction has particle- 
hole symmetry, while the W > 1 states do not. This particle- 
hole symmetry C can also be used as a symmetry protection 
via the Z 2 x Z 2 group generated by {77, C} or by the time- 
reversing symmetry UCt. This fact is well known in the con¬ 
text of the AKLT state, where 77 and C are represented in 
the spin-language as tt rotations about the 2 : and x axes, and 
IICt is the time-reversing symmetry iSyK that flips all com¬ 
ponents of the spins. 

In the context of the argument laid out in the conclusion, 
it seems that particle-hole symmetry can play the same role 
as inversion symmetry in ensuring the edge remains at half¬ 
filling. 


Appendix C: Variants on the HFBI wavefunction 

1. Tuning soft-core bosons to hard-core 

In Equations (3) and (4), the tensor D can be replaced by a 
more general form 


D 


P,*0*1*2 


dp : p = io + ii + i 2 

0 : else 


(Cl) 


which the coefficients dp = 1, 1, v/2, v/6 forp = 0, 1, 2, 3 
in the soft-core state and dp = 1, 1, 0, 0 for p = 0, 1, 2, 3 
in the hard-core state. We can continously tune the coeffi¬ 
cients d 2 and d^ from the soft-core to the hard-core values. 
Upon doing so, we find that the transfer matrix spectrum re¬ 
mains gapped, with the correlation length monotonically in¬ 
creasing from the soft-core state to the hard-core state. Fur¬ 
thermore, the low energy parts of the entanglement spectrum 
do not change significantly through this tuning. Therefore we 
expect that the hard-core and soft-core phases can be adiabati- 
cally connected with a path of local Hamiltonians, and all SPT 
results that apply to one state apply to the other. By choosing 
appropriate values of d 2 and ds, we can also make a state that 
is equivalent to replacing the vacuum 10) in Equation ( 1 ) with 
a constant background of N bosons on each site, N ^ 00 , 
and applying boson annihilation instead of creation operators. 
We can also make a state of spin-5' spins, which is however 
not SU(2)-invariant, where Equation (1) becomes 


IV')=n (E ^*^) n I ^*"=-‘^) • (C2) 

O \*gO / * 

Here, the hard-core state would most closely correspond to a 
state of 5 = 1/2 spins, while the soft-core state corresponds 
to a state of 5 = 3/2 spins. All of these states have the same 
symmetry protection properties. 


2. Interpretation as a Fermionic Wavefunction 

We can also interpret the hard-core variant of the HFBI as 
a wavefunction for spinful fermions on the honeycomb lattice 
at half filling. Note that including the spin, ‘full filling’ of 
a site corresponds to a pair of fermions on each site, so half 
filling occurs with exactly one fermion per site, corresponding 
to two fermions per unit cell. Assuming no spin polarization, 
there must be an equal number of ‘up’ and ‘down’ spins. We 
can bind pairs of opposite-spin fermions into a Cooper pair, 
which yields one Cooper pair per unit cell. As a Cooper pair 
is equivalent to a hard-core boson, we may place the Cooper 
pairs into the hard-core variant of the HFBI. This is equivalent 
to the wavefunction 

i'i^^)=nfE"Ei) io)- (C 3 ) 

o ViGO / 

As the Cooper pair is in a spin singlet state, this wavefunc¬ 
tion preserves SU (2) spin symmetry, in addition to the lattice 
and 7/(1) charge conservation symmetries. It is therefore a 
symmetry-preserving wavefunction of spinful fermions (i.e., 
electrons) on the honeycomb lattice at half filling. However, 
it is not a valid wavefunction for a pure 57/(2) symmetric 
spin model, as it has a vanishing projection onto the subspace 
where each site has exactly unit occupancy. Note that the ne¬ 
cessity to have ‘preformed pairs’ that can then be put into a 
hard-core boson state vividly illustrates the fundamentally in¬ 
teracting nature of this fermionic state. 


3. Inversion Protected Phase 


Additionally, the tensor W in Equation (2) can be replaced 
by the more general form 


W 


ni ...**6 _ 


Px • ^X — 1; — 0 V 2/ 7^ ^ 


I 0 : else 

which corresponds to modifying Equation ( 1 ) to 


(C4) 


li’e) = n E^*^* 1^^- 

O ViGO / 


(C5) 


This does not in general preserve the rotational symmetry 
of the state, but it does if the coefficients Po, ■ ■ ■P 5 are in an 
angular momentum mode 


where £ G (0, 27r/6,..., 5 • 27r/6}. These 6 discrete solutions 
can’t be continously tuned to one another while preserving all 
the lattice symmetries. 

The state |'0^=7r) can be shown to be related to state 
discussed in the main text by a on-site unitary operator 7/ (tt), 
where 

[/(p) = (C6) 

jGB 
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Group 

Generators 

Invariant 

i 



VxVi = 

-I 

— 


{^y} 

Vx Vi = 

■^V ■‘-V 

-I 

- 

Z2 X Zf ^ 

{n, tHT} 

VnVxnxVh^V-yj, = 

-I 


Z2 X Zf 

{n,rnxy} 

VnVrnXyVn^V^yXy = 

-I 


Z2 X Zf ^ 

{nT,,rnT} 

VnXxVxnxV-lV-y^ = 

-I 


Z2 X Zf 

{nix, Tilly} 

VnXxVxnXyVnlV-y^^ = 

-I 



Table II. Summary of symmetry protecting invariants found for the 
state. The degenerate entanglement spectrum cannot be split 
unless all 6 of the minimal protecting symmetry groups are broken. 


Due to this relation, |t/’r= 7 r) and \'4>e=o) have identical correla¬ 
tion lengths and entanglement spectra. However, the protect¬ 
ing symmetries from Table I are mapped using conjugation 
by [/(tt) into a new set of protecting symmetries, shown in 
Table II. Notably, since 

[/(7r)i7IC/(7r)1' = I, (C7) 


this state has doubly degenerate entanglement spectra on odd 
cylinder sizes protected by lattice inversion symmetry alone. 
Thus while the entanglement degeneracy in the HFBI state 
|■^/'r=o) is not split under a staggered field 


H' = h^J2{-iy {b, + bl) 

i 


with (-1)* 


1 i G A 
-1 iGB 


(C8) 


(which fully breaks (7(1) charge symmetry and inversion but 
not the combined symmetry III), the entanglement degener¬ 
acy in the state \'>pi=Tr) would be unsplit by a uniform field, 
which may be physically more interesting. 

A similar mapping for 1-D inversion protected states is 
discussed in Appendix A of Ref. 45. As discussed in Ap¬ 
pendix B, the state \4>t=o) on the W = 1 cylinder is adiabat- 
ically connected to the 1-D Haldane insulator state [21, 45]. 
Correspondingly, the state \4’t=TT) on the W = 1 cylinder is 
adiabatically connected to the 1-D AKLT state. 
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